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ie An abbreviated version of this protocol was published in eLIFE in Nov 2021 
Single-nucleus transcriptomic analysis of human dorsal root ganglion neurons 


DOI: 10.7554/eLife.71752 


Detailed protocol 


Isolation of neuronal nuclei from human DRGs (or tough to dissociate tissues): 


e Freeze tissues on dry ice and store at -80°C until use. Alternatively, cut tissues into small pieces (1-2 mm) and incubate in RNAlater (ThermoFisher, 
cat# AM7021) overnight at RT 
e Remove excess RNAlater and freeze tissue pieces in a microfuge tube on dry ice. Store at -80°C until use 
© Onday of 10x run: prepare homogenization buffer (250 mM sucrose, 25 mM KCI, 5 mM MgCl2, 10 mM Tris, pH 8.0, 1 uM DTT, 0.1% Triton X-100 
(viv). Cool onice. Stock solutions can be prepared ahead of time. 
* Cool down all necessary equipment and solutions. Cool Spectrum™ Bessman Tissue Pulverizer (ThermoFisher, cat#08-41 8-3) in liquid nitrogen or 
dry ice 
e Add frozen tissues to Spectrum™ Bessman Tissue Pulverizer (ThermoFisher, cat#08-4 18-3) in liquid nitrogen or dry ice 
e Pulverize tissues and transfer to glass dounce homogenizer (Fisher Scientific, cat#357538) containing 1ml of cold homogenization buffer 
e Homogenize tissues in glass dounce homogenizer in 1ml of cold homogenization buffer: 5 strokes with “loose” pestle and ~10-20 strokes with “tight” 
pestle ... onice 
e Filter through 40 um cell strainer (ThermoFisher Scientific, cat# 08-771-1) 
e Transfer to microfuge tube (low bind; Sorenson BioScience, cat# 11700) and spin @ ~300-800g at 4° for 5-8 mins (determine what speed and time 
are best for your tissues) 
e Remove supernatant and resuspend pellet in 500 pl of PBS + 1% BSA + SUPERaseln RNase Inhibitor ( 0.2 U/u; ThermoFisher Scientific, Cat# 
AM2696) 
e Incubate on ice for 10-15 mins 
e Add rabbit polyclonal anti-NeuN antibody (Millipore, cat#ABN78) 1:4000 to 1:5000 
e Incubate with rotation at 4° for 30 mins 
e Spin @ 300-800g at 4° for 5-8 mins 
e Remove supernatant and “wash” by adding 1 ml PBS + 1% BSA + SUPERaseIN RNase inhibitor 
© Spin @ ~300-800g at 4° for 5-8 mins 
e Resuspend pellet in 80 pl of PBS + 0.5% BSA + 2mM EDTA. Volume is for 10’ total cells. Adjust volumes accordingly if you have more total cells. 
e Add 20 ul of anti-rabbit IgG Microbeads (Miltenyi Biotech, cat# 130-048-602). Volume is for 10’ total cells. Adjust volume of microbeads accordingly if 
you have more total cells. 
e Incubate at 4° 15-20 mins 
e “wash” by adding 1 ml of PBS + 0.5% BSA + 2mM EDTA 
© Spin @ ~300-800g at 4° for 5-8 mins 
e Remove supernatant and resuspend in 0.5-1 ml PBS + 0.5% BSA + 2mM EDTA 
e Load onto LS column (Miltenyi Biotec, cat# 130-042-401) — follow instruction from manufacturer: 
© put column on magnetic MidiMACS separator (Miltenyi Biotec, cat#130-042-302; MACS MultiStand (cat#130-042-303), add ice cold 3 ml of 
buffer (PBS + 0.5% BSA + 2mM EDTA) to equilibrate 
© Add nuclei to column 
° Wash column with 3 ml ice cold buffer (3x) 
© Remove column from magnet and elute in 5 ml of same buffer using plunger that comes with the LS column. Elute on ice into 15 ml Falcon tube. 


e Spin @ 500g at 4° for 10 mins 

e Remove supernatant 

e Resuspend in 1.5 ml of PBS + 1% BSA + SUPERaseln RNase Inhibitor. 

© Spin @ 300-800g at 4° for 5-8 mins 

e Resuspend in 1-1.5 ml of PBS + 1% BSA + SUPERaseln RNase Inhibitor 

e Strain through 35 um cell strainer (Falcon 352235) to remove any nuclei clumps. 
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Alternatively, use Ultra-Turrax (Laboratory Supply Network, Inc., cat#IKA:3737001) on setting 1 for 10-45 sec on ice (adjust time accordingly ... depending on 
how many nuclei you have) 


o 1% count: Stain 5-10 ul with trypan blue and count (using a hemocytometer) also check for clumps 


e Transfer to low bind microfuge tube and spin @ ~300-800g at 4° for 5-8 mins 

e Remove supernatant and resuspend in desire volume with PBS + 1% BSA + SUPERaseln RNase Inhibitor ... volume based on the first count. 

e Usually | try to go for 1000 nuclei/ul (if there are enough nuclei). | also assume that the first count is not that accurate (too diluted), so usually 
resuspend nuclei in % the calculated volume. 

e Stain 5-10ul with trypan blue and doa 24 count. Adjust volume to desire nuclei concentration. 

¢ Optional: doa 3" count. This should be close to the 2™ count. 

e Load appropriate number of nuclei onto 10X chip. 


Homogenization Buffer: 
250 mM sucrose, 25 mM KCI, 5 mM MgCI2, 10 mM Tris, pH 8.0, 1 uM DTT, 0.1% Triton X-100 


** Can also use EZ PREP buffer (Sigma, Cat #NUC-101) as the homogenization buffer. 


This is a modification of the protocol described in Nguyen et al (2019) Stereotyped transcriptomic transformation of somatosensory neurons in response to 
injury eLife 8:e49679 


Data analysis 
10x chromium data were mapped using CellRanger to a pre-mRNA modified human genome (GRCh38.v25.premRNA). Data analysis used the Seurat V3 


packages developed by the Satija lab and followed standard procedures for co-clustering (Stuart et al., 2019). For sn-RNA sequencing experiments cell 
filtering was performed as follows: outliers were identified and removed based on the number of expressed genes (500—10,000 retained) and mitochondrial 
proportion (<10% retained). Normalization and variance stabilization used regularized negative binomial regression (sctransform). After initial co-clustering 
of data from the different preparations, non-neuronal cell clusters were identified by their gene expression profiles. Clusters not expressing high levels of 
neuronal or somatosensory genes like SVAP25, SCN9A, SCN10A, PIEZO2, NEFH, etc. but instead expressing elevated levels of markers of non-neuronal 
cells including PRP1, MBP, QKI, LPAR1, and APOE were tagged as non-neuronal and were removed to allow reclustering of ‘purified’ human DRG 
neurons. A total of 1837 human DRG neuronal nuclei were included in the analysis (Table 1). The mean number of genes detected per nucleus was 2839 
(range 501-9652), with a standard deviation of 1917. Doublet detection was performed on the individual datasets using DoubletFinder (McGinnis et al. 
2019). For the clustering shown in the main figures the small number of potential doublets (Figure 1—figure supplement 2E) were not removed; principal 
components (PCs) were determined from integrated assay data and PCs 1—16 were used both for UMAP display of the data and for determining clusters. 


The resolution for clustering used ais low stringency (2.0) and closely related clusters without distinguishing markers were merged. The plots shown in 
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